Rapid convergence of time-averaged frequency in phase synchronized systems 
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Numerical and experimental evidence is presented to show that many phase synchronized systems 
of non-identical chaotic oscillators, where the chaotic state is reached through a period-doubling 
cascade, show rapid convergence of the time-averaged frequency. The speed of convergence toward 
the natural frequency scales as the inverse of the measurement period. The results also suggest an 
explanation for why such chaotic oscillators can be phase synchronized. 
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I. INTRODUCTION 

The rich collective behavior, including mutual entrain- 
ment and self-synchronization, in systems of coupled os- 
cillators has been considered by several investigators in 
the past few years (see, for example, Refs. 0, 0, 01 an d 
references therein). Recently, a considerable amount of 
research has been devoted to the study of coupled chaotic 
oscillators and, in particular, to the phenomenon of phase 
synchronization. Provided that the phase can be defined 
fj, E|, two coupled nonidentical chaotic oscillators are 
said to be phase synchronized if their frequencies are 
locked but their amplitudes are not Q, |(| . This appears 
to be a general phenomenon and it has been observed in 
systems as diverse as electrically coupled neurons 0, @ , 
biomedical systems @, chemical systems and spa- 
tially extended ecological systems ^]]- Moreover, the 
potential role of phase synchronization in brain functions 
has been explored [13, 03 ■ 

The most common theoretical approach to phase syn- 
chronization of chaotic oscillators is based on an analogy 
with the evolution of the phase of a periodic oscillator in 
the presence of external noise . This approach leads to 
the conclusion that the dynamics of the phase is generally 
diffusive and the phase performs a random walk. How- 
ever, the effective "noise" in such a description cannot be 
considered as a Gaussian 5-correlated noise in all circum- 
stances. In this paper we show that the effective "noise" 
exhibits strong temporal (anti-)correlations for a gen- 
eral class of chaotic attractors. In particular, we present 
evidence from simulations and experiments that many 
phase synchronized systems of chaotic oscillators where 
the chaotic state is reached through a period-doubling 
cascade show a rapid convergence of the time-averaged 
frequency. The speed of convergence toward the natu- 



ral frequency scales as the inverse of the measurement 
period. This implies that short measurement times may 
suffice for reliable determination of frequencies in those 
systems. 

The outline of the paper is as follows: Section [H] 
presents a sketch of the theoretical background of phase 
dynamics in chaotic systems and the expected scaling of 
the time-averaged frequency with the observation time. 
The results of numerical simulations of globally coupled 
arrays of Rossler oscillators are presented in Sec. 11111 
where features of the chaotic attractor leading to rapid 
convergence of the time-averaged frequency are identi- 
fied. Section llVl contains a description and analysis of ex- 
perimental results on a globally coupled array of electro- 
chemical oscillators where rapid frequency convergence is 
observed. 



II. THEORETICAL BACKGROUND 

An appropriate definition of the phase for chaotic self- 
sustained oscillators can be obtained from the Poincare 
map of the flow. Such a map can be constructed if a 
surface of section in the phase space of the autonomous 
continuous-time chaotic system exists which is crossed 
trans versally by all trajectories of the chaotic attractor. 
Then, for each piece of a trajectory between two crossings 
of this surface, we define the phase as a linear function 
of time: 



(j>(t) = 2tt 



t tn 
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for t n <t < tn-fi- Here, t n is the time of the nth crossing 
of the surface of section. The definition is ambiguous 
because it depends on the choice of the Poincare surface. 
Yet, any choice of a phase variable for chaotic oscillators 
investigated in this paper leads to the same macroscopic 
behavior 

With this phase definition, the phase dynamics can be 
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described by 



A n+1 = M(A n ), (2) 
d(j)/dt = uj(A n ) = L} + r){A n ), (3) 

where the amplitude A n is the set of coordinates of the 
phase point on the Poincare surface at the nth intersec- 
tion and M. defines the Poincare map that takes A n to 
A n+ i. The "instantaneous" frequency u>{A n ) — 2it/T n is 
determined by the Poincare return time T n = t n +i — t n 
and depends in general on the amplitude. Assuming 
chaotic behavior of the amplitudes, we can consider the 
term ui{A n ) to be the sum of the average (natural) fre- 
quency loq and some effective "noise" rj{A n ) with zero 
mean, although this stochastic term has a purely deter- 
ministic origin Thus, Eq.Q has the solution 



<t>{t) = + u Q t + / rj(T)d,T 



(4) 



The variance of the integral J * r\ (t)cZt is given by 

^j\(r)dT^ \ = 2 J\r\t-T')K{T') 

w 2t / dr'K(T') ee 2tD v . (5) 
Jo 

Here, K(r') = (??(t)t?(t 4- r')) and the average (• • • ) is 
taken over the invariant measure of the chaotic attractor. 
The approximate equality holds provided that r/(r) has a 
finite correlation time t c such that K(t') » for r' > t c 
and one considers times t^>t c . This shows that the dy- 
namics of the phase is diffusive and the phase performs 
a random walk as long as t c is finite and D n 0. How- 
ever, as pointed out by Pikovsky and co-workers Q, the 
dynamics of the phase generally differs from stochastic 
Brownian motion because the effective "noise" cannot be 
considered as a Gaussian white noise process. This was 
recently confirmed for a system of locally coupled Rossler 
oscillators where temporal correlations were shown to ex- 
ist ^i). These correlations determine the speed of con- 
vergence of the time average u>(T) = T~ x (4>(T) — <j)(Q)) to- 
ward ojq. The speed of convergence — as measured by the 
standard deviation of the ensemble distribution of u>(T) 
— scales as 1/T. This is in contrast to what one would 
expect if rj were a white noise process where one obtains 
l/VT scaling since uj (T) - luq = 1/T / dtr](t) and the 

standard deviation of J Q T dtrj(t) scales with \/T. Here, 
we provide further numerical and experimental evidence 
that such a scaling generally occurs for oscillators where 
the chaotic state is reached through a period-doubling 
cascade. In particular, the type of coupling between the 
chaotic oscillators seems to be irrelevant. Moreover, we 
show that the temporal correlations induce an alternat- 
ing behavior of K(t'). This leads to an extremely small 
value of D v for these systems and suggests that they can 
generally be phase synchronized. 




T/T„ 



Figure 1: Standard deviation <j^(T) of the u>(T) distribution 
in the phase synchronized state describing the speed of con- 
vergence toward luq. The data points for K — 0.2 are shifted 
down by one decade for clarity. The Poincare surface was 
chosen as the surface spanned by the negative xi axis and the 
X3 axis. The solid lines are best fits to the respective data 
and decay as 1/T a . Note that To = 2-k/ojq depends slightly 
on K. 



III. GLOBALLY COUPLED ROSSLER 
OSCILLATORS 

We consider a L x L array of globally coupled, non- 
identical chaotic Rossler oscillators: 



dx(r,t) 
dt 



= R(x(r,t))- 



-X]T(x(f,i)-x(r,i)), (6) 



where R\ — —uo(r)x2 — £3, R2 — u{v)xi + 0.2x2, R3 = 
X1X3 — 5.9x3 + 0.2. The sites of the lattice are labelled by 
r, K = K' /L 2 , K' is the coupling constant and N is the 
set of all sites on the lattice. We take L — 64 and choose 
the w(r)'s randomly from a uniform distribution on the 
interval [0.99, 1.01]. This ensures that all oscillators are 
chaotic. We find that the speed of convergence of u>(T) 
toward loq in the phase synchronized state scales as 1/T 
(see Fig.^l. This is also true for higher u)(r) dispersion. 

The 1/T scaling observed in the globally coupled 
Rossler system follows from the fact that the local chaotic 
attractors, although one-banded in the phase synchro- 
nized state (see Fig- EJ , have internal structure. 

In the absence of coupling (K — 0) the local non- 
identical Rossler oscillators display a variety of banded 
chaotic attractors because of period-doubling cascades 
and period-3 windows in their vicinity. In the phase syn- 
chronized state for K = 0.06, the local chaotic attrac- 
tors are very similar to one another and have an internal 
structure similar to that of a merged two-banded attrac- 
tor. Although this structure is not obvious in the left 
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Figure 2: Left panel: xi — X2 projection of the local attractor 
at a single site of the lattice for K = 0.06 in the phase syn- 
chronized state. All local attractors look similar to the one 
shown here. Right panel: A typical To = 2-k/ljq histogram 
h(To) of the nonidentical chaotic oscillators for K = 0, i.e., 
the distribution of the natural periods. The exact shape of 
the distribution depends on the w(r) realization. 




Figure 3: Superposition of eight local next amplitude maps 
for K = 0.06. The amplitude A n< \ chosen here is the nega- 
tive x\ coordinate of A n because the X3 coordinate is almost 
constant and very close to zero. The dashed lines correspond 
to the intersection point of the bisectrix with the map. Inset: 
Density distribution p(l) along the curve formed by the next 
amplitude maps. 



panel of Fig. [21 it can be seen clearly in the local next 
amplitude maps (Fig. [3J| and in the density distribution 
p(l) (inset in Fig. ■ The intersection point of the map 
with the bisectrix defines four regions in the map such 
that the two-bandedncss of the attractor can be charac- 
terized by the percentage p of points in the upper left 
and lower right quarter-planes. Indeed, only very few 
points (1 — p = 3%) — typically those for small A n ys 
which are in the lower left quarter-plane of the amplitude 
map — stay in the same region of the bimodal density 
distribution, implying that the majority of the iterates 
map regions of high density on the left to regions of high 
density on the right and vice versa. This nearly two- 
banded character is also seen in the series of T n 's (see 
Fig. 0} . An alternation of long and short return times 
can be observed with rare exceptions. This shows the 
close relation between topological characteristics (next 
amplitude map) and temporal evolution (T„,T„ + i map) 
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Figure 4: Time series of T n for K — 0.06 (left panel) and 
the corresponding T n ,T n +i map (right panel). The arrow 
highlights a deviation from the alternating behavior. These 
deviations occur with 1 — p = 2% at this particular site of the 
lattice which is the percentage of points in the upper right 
and lower left quarter-plane in the lower panel. Note that 
only a short segment of the T n series is shown. 



and suggests the following explanation for the observed 
scaling: 

Assuming that each internal "band" of the attractor 
is associated with a distinct mean frequency wi or u>2 
such that ujq = 2cj]W2/(&>i + ^2)1 an d the system alter- 
nates between bands, the "random" force in the evolution 
equation © for the phase of a single oscillator takes the 
form 



m = (-1) 



(wi + w 2 ) 



(7) 



where k = 1(2) when the system is on band 1(2). The 
noise term fj(t) accounts for deviations of the period 
within a band and deviations from the strict alternations 
between bands. If T is the time needed for n oscillations, 
it follows from Eq.Q that 



w(T) - UJq = 



n even. 



(8) 



where K = 1(2) corresponds to band 1(2) appearing first 



in the series. The term 



2tt(u 



scales with 1/T, while 



(wi+w 2 )T 

the term J Q dt?4p- scales with l/\/T provided that fj(t) 

is (^-correlated noise. Since the amplitude of the 1/vT 
term is very small, two scaling regimes can be identi- 
fied for the maximal deviation from luq: A 1/T scaling 
on intermediate time scales which is dominated by the 
switching between the "bands" and a 1/y/T scaling on 
long time scales. The intermediate time scale scaling be- 
havior can indeed be observed for each oscillator: The 
envelope of Q(T) — loq decays as 1/T. Note that the 
standard deviation of the ensemble distribution of Cj(T) 
decays as 1/VT for identical chaotic Rossler oscillators 
and K — — as expected from the law of large numbers. 
In the case of nonidentical oscillators, phase synchroniza- 
tion implies that the phases are synchronized such that 
yf ((u>(T) — w ) 2 ), where the average is taken over the 
ensemble, adopts the scaling of a single oscillator repro- 
ducing the observed scaling and explaining the results. 
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Figure 5: Superposition of the next amplitude maps for K — 
0.2 at the same sites as in Fig. [3] Inset: Normalized point 
density along the curve formed by the next amplitude maps. 
Three regions of high density can be identified. 



This result does not depend on the type of coupling of 
the oscillators as confirmed by the numerical results pre- 
sented here and in Ref. 0|. This explanation is 
based exclusively on the fact that the series of the in- 
stantaneous frequencies alternates with the exception of 
rare events and that the system is phase synchronized. 

The above argument easily generalizes to attractors 
with multiple internal "bands" which we observe in the 
phase synchronized state for higher coupling. Figure 
shows that the deviations from a three-banded structure 
are rather small for K = 0.2. Three bands A, B and C 
can be identified such that A is mapped to B and C to 
A. However, band B is mapped to C and Z, and Z back 
to B such that a pure three-banded structure does not 
exist. The invariant measure of the map in area Z is very 
low (approximately 10%). 

The alternation of positive and negative values for rj{t) 
implies that K(t') alternates as well. Hence, D v is very 
small, i.e., a high degree of phase coherence persists. Our 
findings suggest that this should be generally expected 
for oscillators where the chaotic state is reached through 
a period-doubling cascade. Thus, these oscillators would 
be phase coherent suggesting that they can generically 
be phase synchronized. 

The change in the structure of the local attractors as 
the coupling strength is changed is the analog of the bi- 
furcation structure seen in homogeneous locally coupled 
Rossler oscillators where one finds period doubling dif- 
ferent from the isolated Rossler oscillator [l6j| . 



IV. GLOBALLY COUPLED 
ELECTROCHEMICAL OSCILLATORS 

We now study the convergence of the time-averaged 
frequency in an array of globally coupled electrochemical 
oscillators. A standard three electrode electrochemical 



cell consisting of a nickel working electrode array (64 1- 
mm diameter electrodes in an 8 x 8 geometry with 2 mm 
spacing), a Hg/Hg 2 S04/K 2 S04 reference electrode and 
a Pt mesh counter electrode were used. The potentials 
of all of the electrodes in the array were held at the same 
value (V = 1.310 V) with a potentiostat (EG&G Prince- 
ton Applied Research, model 273). Experiments were 
carried out in 4.5 M H2SO4 solution at a temperature of 
11 °C. The working electrodes were embedded in epoxy 
and reaction takes place only at the ends. The currents 
of the electrodes were measured independently at a sam- 
pling rate of 100 Hz and thus the rate of reaction as a 
function of position and time was obtained. 

The electrodes were connected to the potentiostat 
through 64 individual parallel resistors (R p ) and through 
one series collective resistor (i? s ). We employed a method 
of altering the strength of global coupling while holding 
all other parameters constant. The total external resis- 
tance was held constant while the fraction dedicated to 
individual currents, as opposed to the total current, was 
varied. A total resistance can be defined as 

Rtot = R s + i? p /64 . (9) 

In these experiments i? to t = 14.2 fi. The series resistor 
couples the electrodes globally. The parameter e, the 
ratio of series to total resistance, is a measure of the 
global coupling 



For e = 0, the external resistance furnishes no additional 
global coupling; for e = 1, maximal external global cou- 
pling is achieved. 

The population of chaotic oscillators is characterized 
by certain amount of heterogeneity [Hj • Without added 
coupling there is a distribution of frequencies of the oscil- 
lators with a mean of 1.219 Hz and a standard deviation 
of 18mHz. With weak added global coupling (e = 0.1) 
a nearly phase synchronized state occurs in which 63 of 
64 oscillators have a frequency of 1.230 Hz, and that of 
the remaining element has a frequency of 1.237 Hz. The 
results are shown for this region of high phase synchrony. 

Experimental Results 

The chaotic attractor of a representative single clement 
of the coupled system is shown in Fig. HJ] The chaotic 
state is adjacent to a period-three window and a period 
doubling sequence. The reconstructed attractor (Fig. EJi) 
is low-dimensional (correlation dimension 2.3 ± 0.1). To 
obtain the phase, the Hilbert transform approach is ap- 
plied 0. The two-dimensional embedding using Hilbert 
transform is shown in Fig. |SJd. As in the case of globally 
coupled Rossler oscillators, it is difficult to distinguish 
the banded character of the attractor. 

The merged banded structure is more clearly seen in 
return maps obtained from the series of return times 
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Figure 6: Experiment. Chaotic dynamics of a single element 
at e = 0.1. a. Reconstructed attractor using time delay 
coordinates, b. Two dimensional phase space reconstruction 
using the Hilbert transform. The solid circle represents the 
origin used for phase calculations. 
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Figure 7: Experiment. Time series and maps of discretized 
chaotic dynamics of a single element at e — 0.1. a. Time series 
of return time (T„) of the oscillations, b. One dimensional 
map using the return time. The dashed lines correspond to 
the average period of oscillations (To = 0.81 s) which in this 
case is identical with the intersection of the map with the 
bisectrix. 



(T n ). The series of the return time clearly shows an al- 
ternation of long and short oscillations (Fig. EJi). The 
map constructed using T n exhibits an approximately one- 
dimensional character (Fig. \7jp). 

The oscillators in the population are slightly different. 
Figure |BK shows a superposition of the ID maps of all 
the oscillators. Clearly, the majority of the phase points 
lie in the upper left and lower right boxes indicating two- 
banded character. The probability of two-banded oscil- 
lations varies from one site to the other (Fig. |Hb) with a 
mean of p me an = 0.85 and a standard deviation of 0.05. 
We note that the two-banded character of the uncou- 
pled system of uncoupled oscillators {e = 0, p mea n = 0.7, 
Figs. |Ht and d) is smaller than that of the phase synchro- 
nized region. Therefore, during the transition to phase 
synchronization with increasing the coupling strength the 
two-bandedness observed in the time series of oscillations 
is more pronounced. 

As expected from the numerical simulations and theo- 
retical considerations, the strong (anti)correlations in the 
phase dynamics have pronounced effects on the speed 
of convergence of frequencies. Figure El shows that the 
standard deviation of frequencies is proportional to 1/T 
at e — 0.1. (The fluctuations around the fitted line are 
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Figure 8: Experiment. Characterization of the two-banded 
character of the population of sixty four oscillators. Top row: 
e — 0.1. Bottom row: e = 0. a. Superimposed one di- 
mensional maps of the sixty-four sites. The dashed lines cor- 
respond to the average period of oscillations (T = 0.81s); 
e = 0.1. b. The distribution of probability (p) of two-banded 
oscillations; e = 0.1. c. Superimposed one dimensional maps 
of the sixty-four sites; e = 0. d. The distribution of probabil- 
ity (jp) of two-banded oscillations; e = 0. 




Figure 9: Standard deviation of the experimental frequency 
distribution as a function of the applied time average (the 
thick line is a linear fit with slope -0.97) for e = 0.1. 



caused by slowing down and speeding up within a cycle.) 
Similar scaling results are obtained for other (larger) cou- 
pling strengths. 



V. CONCLUSION 

For many systems of coupled oscillators, the chaotic 
state is reached through a period-doubling cascade. This 
implies that a Poincare surface exists such that the next- 
amplitude map is similar to that in Figs. or [S] As long 
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as the correspondence between the next-amplitude map 
and the relation of subsequent return times holds, this 
should lead to a small value of D v and, in particular, to 
a 1/T scaling in the phase synchronized state. The ex- 
perimental results show that such a structure of the series 
of return times can exists even if no well-behaved next- 
amplitude map can be identified. Our findings further 
suggest that chaotic oscillators where the chaotic state 
is reached through a period-doubling cascade can gener- 
ally be phase synchronized due to the small value of D n , 
i.e., high degree of phase coherence. This is in accord 
with findings in Ref. where different chaotic systems 
were analyzed with respect to their ability to be phase 
synchronized. There it was found that the Lorenz sys- 
tem (where the chaotic attractor is not reached through a 
period-doubling cascade) cannot be phase synchronized. 
In particular, the next amplitude map has a very different 



structure, i.e., a logarithmic singularity. 

Finally, in many biological and physical systems only 
short time series are available. Our results show that 
the determination of frequencies in weakly coupled oscil- 
latory systems may be obtainable from such short time 
series. 
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